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ABSTRACT 

The BL Lacertae (BL Lac) object lES 2344+514 (lES 2344), at a redshift 
of 0.044, was discovered as a source of very high energy (VHE) gamma rays by 
the Whipple Collaboration in 1995 (Catanese et al. 1998). This detection was 
recently confirmed by the HEGRA Collaboration (Tluczykont et al. 2003). As 
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is typical for high-frequency peaked blazars, the VHE gamma-ray emission is 
highly variable. On the night of 20 December, 1995, a gamma-ray flare of 5.3- 
sigma significance was detected, the brightest outburst from this object to-date. 
The emission region is compatible with a point source. The spectrum between 
0.8 TeV and 12.6 TeV can be described by a power law 

,f!^,, = (5.1 ± l.Ost ± 1.2,J X io-7(E/TeV)-2-^"±°-^^-±'^-°^- ^ . 

dE dA dt ^ St sy> \ I ) rpgy ^2 s 

Comparing the spectral index with that of the other five confirmed TeV blazars, 
the spectrum of lES 2344 is similar to lES 1959-1-650, located at almost the same 
distance. The spectrum of lES 2344 is steeper than the brightest flare spectra 
of Markarian 421 (Mrk 421) and Markarian 501 (Mrk 501), both located at a 
distance about 2/3 that of lES 2344, and harder than the spectra of PKS 2155- 
304 and H 1426-1-428, which are located almost three times as far. This trend 
is consistent with attenuation caused by the infrared extragalactic background 
radiation. 

Subject headings: BL Lacertae objects: individual lES 2344-1-514) — gamma 
rays: observations 

1. Introduction 

To date, the only confirmed extragalactic gamma-ray sources at energies > 100 GeV 
(very high energy, VHE) are BL Lacertae (BL Lac) objects and the giant radio galaxy M87 
(Beihcke et al. 2004). BL Lac objects are active galactic nuclei (AGN) with (1) characteristic 
radio/optical/X-ray flux, (2) the absence of emission lines with observed equivalent width 
greater than 5 A, and (3) a CA II "break strength" smaller than 25% (Perlman et al. 1996). 
These criteria define an object with strong nonthermal emission which almost completely 
masks the thermal emission from the surrounding host galaxy. The spectrum, in a vF^, 
representation, shows a double-peaked structure. The only type of BL Lac objects detected 
so far to emit VHE emission are high-frequency peaked. For these objects, the low-energy 
component peaks in the soft to hard X-ray regime and the high-energy component peaks 
in the VHE regime. The six confirmed VHE BL Lac objects are: Mrk 421 (Punch et al. 
1992; Petry et al. 1996), Mrk 501 (Quinn et al. 1996; Bradbury et al. 1997), lES 2344+514 
(Catancse et al. 1998; Tluczykont ct al. 2003), lES 1959+650 (Nishiyama et al. 1999; Holder 
et al. 2003; Aharonian et al. 2003), PKS 2155-304 (Chadwick et al. 1999; Aharonian et al. 
2005), and H 1426+428 (Horan et al. 2002; Aharonian et al. 2002a). The emission level 
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around the two peaks is highly variable, and changes in the spectral shape with flux level 

have been measured for Mrk 421 (Krennrich et al. 2002; Aharonian et al. 2002b; Krennrich 
et al. 2003), Mrk 501 (Djannati-Atai et al. 1999; Aharonian et al. 2001), and lES 1959+650 
(Aharonian et al. 2003). For the other three BL Lacs, variations of the spectral shape with 
flux level have neither been established nor ruled out. 

The VHE observations reported here were carried out by the VERITAS (previously 
Whipple Gamma Ray) collaboration using an imaging atmospheric Cherenkov telescope 
(Weekes et al. 1989). The telescope, of 10 m diameter, is located on Mt. Hopkins at an 
altitude of 2320 m above sea level. At the time of observations, the imaging camera consisted 
of 109 photomultiplier tubes, each viewing 0.259° of the sky and arranged in a closed-packed 
hexagonal pattern. The telescope and the data acquisition are described in Cawley et al. 
(1990). 

The organization of this paper is as follows: The status of observations on lES 2344 is 
summarized in Section 2. The VHE data and analysis techniques are presented in Section 3. 
This is followed in Section 4 by a description of the gamma-ray simulations necessary for 
the spectral reconstruction, including estimation of the gamma-ray energy in Sect. 4.1. The 
measured VHE spectra are presented in Section 5 and are briefly discussed and summarized 
in Section 6. 



2. Observational Status of lES 2344-|-514 

The BL Lac object, lES 2344-^514, at a redshift of 0.044, was detected in the Einstein 
Slew Survey (Elvis et al. 1992) in the energy range 0.2-4 keV. The survey was constructed 
from data collected during the HEAO-2 mission from 1978-1981. lES 2344 was identifled 
as a BL Lac object in Perlman et al. (1996). The non-contemporaneous spectral energy 
distribution of lES 2344 is shown in Figure 1. Observations at all wavelengths show lES 2344 
to be an unresolved point source. The central black-hole mass is lO^-^'^='='^-^^ M©, derived from 
stellar velocity dispersion measurements (Barth et al. 2003). In the optical regime, a point 
source with an underlying elliptical host galaxy can be fitted with a radius (half-width at 
half-maximum) of re = 7.12 ±0.02 kpc {Ho=50 km s"^ Mpc"^ and qo=0) (Urry et al. 2000). 

The optical and far-infrared emission from lES 2344 contains significant contributions 
from the host galaxy. The total photometry by the 2 Micron All Sky Survey (Jarrett et al. 
2003) and by the Hubble Space Telescope (HST) (Urry et al. 2000), labeled "Galaxy light" 
in Figure 1, lie well above the value expected from pure synchrotron emission in the jet. Ob- 
servations with the HST in 1996 measured a R-band brightness of the nucleus of 16.83±0.05 
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mag from a fit of a point source plus galaxy convolved with the point spread function of the 
telescope (Urry et al. 2000). During continued monitoring through 1998, the R-band bright- 
ness varied between 16.47 mag (Nilsson et al. 1999) and 17.00 mag (Falomo and Kotilainen 
1999), indicating optical variability. An optical monitoring program in 2000/1 by Xie et al. 
(2002) found short time scale variability to be weak, with maximum intraday variability of 
AV = 0.18 mag, AR = 0.1 including galaxy light. A relatively large brightness decrease of 
0.35 mag was observed in the V-band over 2 weeks in January 2001. 

lES 2344 showed X-ray variability on the time scale of a few hours in the 0.1 - 10 keV 
energy band during a week-long campaign in 1996 using the BeppoSAX satellite (Giommi 
et al. 2000). A follow-up observation in 1998 found lES 2344 to be in a very low state, 
implying a frequency shift by a factor of 30 or more of the peak synchrotron emission. 
They suggested the interpretation that two distinct electron populations contribute to the 
synchrotron emission; one steady low-energy component, the other producing soft to hard 
X-rays with rapid time variability. 

lES 2344 has been monitored by the Whipple Collaboration since 1995 (Catanese et al. 
1998). Recently, the HEGRA Collaboration reported an independent confirmation of this 
source (Tluczykont et al. 2003). On the night of 5 December, 1996, Whipple VHE and 
BeppoSAX X-ray observations overlapped for 28 minutes, for which wc show the X-ray 
spectrum and VHE flux upper limit in Figure 1. The 99.9% VHE flux upper limit at 
energies > 350 GeV was calculated as in Catanese et al. (1998). 

In the VHE band, the object was observed in a flaring state during the night of 20 

December, 1995, with a significance of 5.3 a, the strongest gamma-ray flare measured from 
this object to date. The quiescent flux level of lES 2344, compared to the flare presented here, 
is about 50 times lower (Tluczykont et al. 2003). The detection of VHE gamma rays from 
lES 2344 in December 1995, reported by the Whipple Collaboration at the 1997 International 
Cosmic Ray Conference (Catanese et al. 1997), was considered tentative because follow- 
up observations by this and other VHE observatories through 1997 did not detect further 
evidence for gamma-ray emission. Monitoring by the Whipple Collaboration from 1998 to 
2000, however, showed again a small positive excess (Badran 2001). A summary of all 
published VHE observations of this source is given in Table 1. An initial measurement of the 
VHE gamma-ray spectrum covering the entire 1995/6 observing season yielded a spectrum 
of (1.14 ±0.50) X io~'''ii;~2-29±o.43 'jgy-i ]^-2 g-i^ statistical error only, over the energy range 
Q.5<E < 5.0 TeV with x^/ndf = 3.2/2 (Bussons- Gordo 1998a,b). 

In Figure 2, we show the two-dimensional gamma-ray sky map during the flare. The 
gamma-ray map was constructed from a partial data set, referred to as 'B' in Section 4. 
The emission region is compatible with a point source. This was determined using a Monte- 
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Carlo simulation of the telescope response to a point source of gamma rays. A small telescope 
pointing error of less than 0.05° may have been present during the observations, but due to 
the lack of bright stars in the field of view, we are not able to determine this in retrospect. 
The centroid of the measured gamma-ray emission is displaced from the known location of 
lES 2344 by RA 0.02 ± 0.02° and DEC 0.03 ± 0.02°. With in a conservatively estimated 
0.1° circle of confusion are located three galaxies and two radio sources, but no other X-ray 
sources. Thus, the gamma-ray emission likely originated from lES 2344. 

The EGRET 95% confidence level upper limit for lES 2344 is 6.98 x 10"^ cts cm^^ s^^ 
E > 100 McV (Hartman ct al. 1999). The peak response for most sources detected with 
EGRET lies at around 300 McV, this corresponds to an upper limit at 300 MeV of about 
3.4 X 10-^^erg cm^^ 

3. Gamma- Ray Flare and Background Data 

Observations with the 10 m telescope were carried out in two pointing modes: (1) with 
the source in the center of the field of view (ON observation) and (2) with the telescope 
pointing offset from the source direction by 30 minutes in RA, called OFF observation. The 
OFF observation is a measurement of the background caused by cosmic rays. On the night 
of 20 December, 1995, four ON observations were taken with a combined exposure time was 
110 minutes. The last ON observation during the night was not complemented with an OFF 
observation, as is necessary for spectral measurements. Therefore, an OFF observation was 
selected from the 1995/6 season based on its similarity to the ON observation in elevation, 
cosmic-ray rate, and night sky brightness. For each ON observations. Table 2 lists the 
UTC start time, the average observing elevation, the throughput factor for both, ON and 
OFF observations, and the measured gamma-ray rate. The throughput factor measures 
the cosmic-ray rate relative to a reference observation taken under clear skies (Lebohec 
and Holder 2003). The weather during these observations was rated "A" by the observers, 
meaning clear skies. Figure 3 confirms this by comparing these observations with other "A" 
weather observations taken between October 1995 and April 1996. 

The standard analysis method for data taken with the 10 m telescope (Reynolds et al. 
1993) includes conditioning of the images, parameterization, and selection of gamma-ray like 
events. Conditioning of the images consists of: (1) fiat- fielding of the relative gain between 
pixels, (2) equalizing the sky brightness between the ON and OFF observation (Cawley 
1993), and (3) removing pixels that arc below a certain signal-to-noise ratio. Images are 
then parameterized by their RMS width and length, their distance from the center of the 
field of view (Hillas 1985), and the orientation angle of their major axis relative to the 
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pointing direction of the telescope, alpha (Weekes et al. 1987). The total amount of hght 
collected is referred to as the size of the image. 

The gamma-ray signal is derived from the excess number of events between ON and 
OFF runs, where only those images are selected that are hkely to have been produced by a 
gamma-ray source located at the center of the field of view. The rate given in Table 2 shows 
the gamma-ray rate after application of one particular set of selection criteria (cuts), called 

Supercutsl995 (Catanese et al. 1998). These cuts are not optimal for spectral measurements 
because the selection efficiency for gamma-rays decreases dramatically with energy. There- 
fore, a different set of cuts was developed empirically using simulated gamma-ray events, 
see Section 4.1. The analysis of the data given in Table 2 is in agreement with a previous 
analysis by Catanese et al. (1998). 

4. Calibration and Spectral Reconstruction 

The data contain a relatively low gamma-ray rate and were taken over a wide range of 
elevations. To obtain an accurate energy calibration our analysis technique requires us to 
analyze different elevation ranges separately. To maintain a good signal-to-noise ratio, the 
data were combined at the two average elevations of 58° and 41° and referred to as datasets 
A and B, respectively. The data in these two sets were taken sequentially during the night, 
allowing us to investigate time variability in the emission level. 

A total of 500,000 gamma-ray initiated showers were simulated at each elevation; these 
were distributed in energy randomly according to a power law of index -2.5 and covering a 
circular area around the telescope axis. The Monte-Carlo simulations of gamma-ray initi- 
ated particle showers in the atmosphere and subsequent detection of Cherenkov photons by 
the telescope were carried out with the Grinnell-ISU (GrISU) package^. At 58° elevation, 
simulations were carried out over the energy range 0.1-100 TeV and impact radius less than 
300 m. At 41° elevation, simulations were carried out over the energy range 0.3-100 TeV 
and impact radius less than 350 m. The low-energy cut-off was chosen to extend beyond 
the range of energies of events that trigger the telescope. The night-sky brightness level of 
simulated showers was matched to that measured from the data. 

The absolute light throughput of the telescope was calibrated with Cherenkov images 
of muons recorded by the telescope. For this, only complete muon rings were selected using 
a specially developed algorithm (Schroedter 2004). This ensures that the total amount of 



^Available at http://www.physics.utah.edu/gammaray/GrISU/ 
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light is well known. The light throughput factor derived in this way was used to measure 

the spectrum of the Crab Nebula during 1995/6 season. With statistical (st) and systematic 
(sy) errors, the spectrum of the Crab Nebula between 0.3 TeV and 13 TeV can be fitted by 
{4.2 ± 0.3 St ±0.7 sy) X io-7^-2-38±o.o8.t±o.o4., rpgy-i ^-2^-1 ^j^^ xLn/ndf = 3.2/(9-2) . This 
is compatible with other measurements (Mohanty et al. 1998; Hillas et al. 1998). 

The energy resolution of the spectral analysis depends on rejecting cosmic-ray images, 
and selecting only those gamma-ray images with well defined image parameters. The col- 
lection area near the triggering threshold is difficult to model in the simulations and hence 
a software cut on the minimum brightness is applied that lies substantially above the hard- 
ware threshold. The following set of loose cuts were then applied to data and simulations: 
0.31°< distance <1.1°, length/size <0.00085 °/dc, max2 > 65 dc, and alpha < 25°. 

The differential trigger rates at 41° elevation, data set B, are shown in Figure 4 for a 
spectrum with differential index of -2.5. The peak trigger rate occurs at an energy of 1.4 TeV 
for spectral cuts, described below, and 2.1 TeV with Supercutsl995. With these cuts, 90% 
of the triggers occur above 1.05 TeV and 1.67 TeV, respectively. The collection area, shown 
in Figure 5, reaches 10% of its maximum value of 170,000 m^ at an energy of about 1.1 TeV 
for spectral cuts. The differential trigger rates at 58° elevation, corresponding to data set A, 
arc shown in Figure 4. The peak trigger rate occurs at an energy of 0.69 TcV for spectral 
cuts and 1.1 TeV with Supercutsl995. For these two sets of cuts, 90% of the triggers occur 
above 0.48 TeV and 0.75 TeV, respectively. The collection area, shown in Figure 5, reaches 
10% of its maximum value of 136,000 m^ at an energy of about 0.51 TeV for spectral cuts. 



4.1. Event Selection and Energy Estimation 

The spectral analysis method has been described in Petry et al. (2002); Mohanty et al. 
(1998). Simulations at 41° elevation show an energy resolution of rms{A\ogE) = 0.15 or 
rms{AE/ E) = 0.40, with an energy estimation bias of |Alogii^| = 0.018 over the energy 
region E = 0.8 TeV to 40 TeV. This energy range begins at 10% of the peak collection area. 
A cut-off at the high energies is necessary as the limited field of view of the camera truncates 
large showers and the estimated and true energies begin to diverge. At 58° elevation, the 
energy resolution is rms{A\ogE) — 0.16 or rms{AE/E) — 0.49) and |Alog£'| = 0.012 over 
the energy range E = 0.4 TeV to 25 TeV. 

The gamma-ray signal is contaminated by a large fraction of cosmic-ray events. To reject 
this background, cuts are imposed on the parameters distance, width., length., and alpha. 
The cuts derived from the Monte-Carlo simulations scale with size so that the efficiency of 
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selecting gamma rays remains unchanged as a function of energy. The fraction of gamma rays 
passing the cuts for simulations at 41° elevation is 86%, and it is 87% at 58° elevation. The 
distributions of the parameters width, length, and alpha are shown for simulated gamma rays 
in Figure 6. The cuts arc chosen at a nominal 2 standard deviations around the mean value. 
The simulations at 58° are limited by statistics at high energies, making the cuts somewhat 
inefficient. In particular, the upturn at large dc value of the alpha-cut is unphysical, but 
the cut-level still remains below the Supercutsl995 value of 15°. The unphysical upturn is 
due to the second order polynomial used in fitting the cut-level. For comparison, the level 
of Supercutsl995 is also shown in Figure 6. 



5. Flare Spectra 

The number of excess gamma-ray events in each energy bin after application of all cuts 
is presented for both data sets in Tables 3 and 4. Due to the very small signal, the bin 
width is chosen at twice the energy resolution A(log£') = 0.3 (Petry et al. 2002). Flux upper 
limits are given if the gamma-ray significance is less than 1 cr in the energy bin. The upper 
limits are at the 98% confidence level and calculated according to the method of Helene 
(1983). The spectra for the two data sets A and B are shown in Figure 7. The error bars 
show the statistical error only. 

For dataset B, the power law fit to the spectrum over the energy range from 0.8 TeV to 
12.6 TeV is given by 

(5.1 ± 1.0,^ ± 1.2,y) X 10-7^-2-54±0.17.*±0.07., ^ ^ 



dEdAdt ' ■ ■ ■ TeVm^s' 

with Xmin/^(^f — -2/(4 — 2). The probability for this data to randomly arise from the 
power-law fit is 0.9. The statistical error represents the 68% confidence interval (CI) for a 
fit with one free parameter while the other parameter is frozen at its optimum value. The 
68% CI with two simultaneous free parameters, defined by xLin + 2. 3, is shown in Figure 8. 

The systematic errors of the flux constant and spectral index arising from the energy 
calibration and the cut-tolerance are indicated in Figure 8 by crosses. The cut tolerance, 
with a nominal value of 2 standard deviations, was varied between 1.5 and 2.5 standard 
deviations; the level of the muon-bascd energy calibration is ±10%. The uncertainty in 
the energy calibration affects mostly the flux constant. For example, a 10% change in the 
energy calibration changes the flux constant by 25% (30%) if the spectrum has a differential 
index of -2.5 (-3.0). In addition, due to the large elevation range covered, a small systematic 
uncertainty on the order of 10-15% is intrinsic to the GrISU simulations (Krennrich et al. 
1999). The spectral index is affected mostly by varying the cut tolerance. It should be 
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noted that the systematic error evaluated in this way is smaller than the statistical error. 
This means that a good estimate of the systematic error is not possible with this method; 
nevertheless it does indicate the relative importance of the two sources of error. 

For data set A, the power law fit over the energy range from 0.4 TeV to 1.6 TeV is given 

by 

= (1.9 ± 0.6,* ± 0.6,,) X io-7£;-3-3±o.7..±o.7., J ^ (3) 



dE dAdt ^ ' ' ' TeV m2 s' 

and the confidence interval contours are shown in Figure 8. 

As the spectral indexes of the two spectra are compatible, it is possible to adjust the flux 
constant of the less significant spectrum (set A) so that it overlaps, in a least-squares sense, 
with the spectrum of set B. However, as the statistical significance of data set A is very small 
compared to set B, combining the two data sets results in an insignificant improvement in 
the statistical error of the spectral index. Therefore, the spectral measurement of lES 2344 
derived here, is best represented by spectrum of data set B, alone. 



6. Discussion 

lES 2344 is a variable source; during the fiare on 20 December, 1995, the gamma- 
ray emission from lES 2344 was about 50 times brighter than during the quiescent phase 
measured several years later. To obtain an accurate energy calibration, our analysis technique 
required us to analyze data taken at different observing elevations separately. Therefore, we 
split the data into two sets, A and B, with 56 and 38 minutes exposure time, respectively. 
The data sets were taken consecutively during the night. The spectral indices measured 
from the two data sets are compatible with each other. The increase of the flux constant 
over the two hours of observation, though also not very significant, is not unexpected as 
very large variability of the VHE fiux on time scale of hours has been observed for other 
blazars (Gaidos et al. 1996; Quinn et al. 1996; Holder et al. 2003; Aharonian et al. 2005). 

The measured VHE spectra are attenuated through pair production with the infrared 
extragalactic background light (EBL) (Nikishov 1962). Due to the EBL spectral shape, 
the attenuation manifests itself as a steepening of the measured VHE spectrum between 
roughly 1 and 5 TeV and becomes more pronounced with larger redshift. A cut-off feature is 
thus expected in the VHE spectra, if variations in the intrinsic VHE spectrum are ignored. 
Such a cut-off feature has been estabhshed for Mrk 421 and Mrk 501 (Aharonian et al. 
2001; Krennrich et al. 2001, 2002). Their spectra can be described with a power law with 
exponential cut-off: dN/dE oc exp (E/Eq). The cut-off energy, £"0, differs between the 
two blazars by 2.6±1.2 TeV (Aharonian et al. 2002b). Unfortunately, for lES 2344 the low 
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st atistical significance of the spectrum precludes the measurement of such a cut-off energy. 

VHE spectra are now available for all six confirmed TeV blazars. The power law spec- 
tral indexes of fits to the brightest flares from the blazars appear to steepen with redshift 
(Schroedter 2005). The spectral index of the lES 2344 VHE flare is steeper than the brightest 
flare spectra of Mrk 421 and Mrk 501, both located at about 2/3 the distance of lES 2344. 
The flare spectra of PKS 2155-304 and H 1426+428, located almost three times as far, 
are softer than that of lES 2344. The spectral index of the lES 2344 flare is similar to 
the flare spectrum of lES 1959+650, which is located at almost the same redshift. This 
trend is consistent with attenuation caused by the infrared extragalactic background radia- 
tion (Schroedter 2005; Stecker 1999). Alternatively, galaxy evolution might be responsible 
for the observed spectral steepening with redshift. For example, if younger galaxies have 
enhanced mid-infrared radiation nearer to the central black hole, then this would produce 
gamma-ray attenuation indistinguishable from that caused by the EBL. 

No contemporaneous measurements at other wavelengths were taken during the gamma- 
ray flare of lES 2344 on 20 December, 1995. This precludes the application of models to 
constrain the gamma-ray production mechanism, because the gamma-ray emission is known 
to be highly variable. Almost one year later, on 5 December, 1996, a simultaneous TeV 
/ X-ray observation occurred together with the BeppoSAX satellite. The detailed X-ray 
spectrum measured during this night (Giommi ct al. 2000) is complemented, however, only 
by an upper limit of the TeV flux, again precluding models to be signiflcantly constrained. 

We acknowledge the technical assistance of K. Harris, T. Lappin, and E. Roache. This 
research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 
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Fig. 1. — The spectral energy distribution of lES 2344 along with the VHE flare spectrum 
obtained with the Whipple 10 m telescope {line with error bars). Also shown is the VHE 
flux upper limit for the night of 5 December, 1996 (open circle). Other data were taken 
from the following sources: 365 MHz from Texas radio survey [filled circle) (Douglas et al. 
1996), 1.4 GHz from Grccnbank (filled circle) (White and Becker 1992), 4.85 GHz from 
Grcenbank [filled circle) (Gregory and Condon 1991), 8.4 GHz from VLA (filled circle) 
(Patnaik et al. 1992), galaxy photometry at millimeter wavelength {filled diamond) (Stevens 
and Gear 1999), galaxy photometry at K, H, and J-bands from 2MASS {line segment) 
(Jarrett et al. 2003), galaxy and nucleus R-band photometry obtained with Hubble Space 
Telescope and corrected for interstellar reddening {filled diamond) (Urry et al. 2000). X-ray 
observation with BcppoSAX {see legend) from (Giommi et al. 2000), upper limit at 300 MeV 
from EGRET {filled triangle) (Hartman et al. 1999). Quiescent VHE gamma ray flux during 
the period 1997-2002 from HEGRA {open diamond) (Tluczykont et al. 2003). 
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Fig. 2. — Gamma-ray sky map of the field of view around lES 2344, with position indicated 
by the red cross. The colors show excess counts with overlaid significance contours in steps 
of one standard deviation per contour. The dotted hues show the RA and DEC. 
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Fig. 3. — Relative cosmic ray rate of lES 2344 flare data compared to observations carried 
out under clear skies (rated "A" by observers) in 1995/6. For clarity, the error bars are only 
shown for the flare data. 
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Fig. 4. — Simulated trigger rate at 41° elevation (left) and 58° (right) of gamma-rays dis- 
tributed with a power law index of -2.5. The two lines show the rate after application of 
spectral cuts (bold) and Supercutsl995 (thin). 
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Fig. 5. — Collection area for gamma rays at 41° elevation {left) and at 58° elevation (right) 
in 1995 for spectral cuts (bold line) and Supercutsl995 cuts (thin line). 
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Fig. 6. — Simulated parameter distributions and cut levels with log{size) after application 
of loose cuts, see text. Top row: simulated at 41° elevation, bottom row. 58° elevation. The 
dots are simulated events and crosses represent the average. The solid lines shows the fitted 
polynomial to the average. Dashed lines show the actual cut chosen at a tolerance of two 
standard deviations around the average. Dotted lines show the cut level of Super cuts 1995. 
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Fig. 7. — Differential flux spectrum of lES 2344 on 20 December 1995. Spectra from data 
sets A (diamonds) and B (circles) are shown together with power law fits {solid lines). 
The shaded regions show the confidence interval of the power law fits and were obtained by 
varying both parameters to their individual 68% confidence interval. 
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Fig. 8. — Confidence regions corresponding to spectrum of data set A {dashed lines) and 
set B (solid lines). Confidence regions are shown with probabihty content of 40%, 68%, and 
90% for the simultaneous values of the spectral index and fiux constant. Also shown are the 
systematic error on the fiux constant and spectral index (crosses) 
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Table 1: Summary of VHE measurements of lES 2344. 



Date 


Reference 


Exposure 




Integral Flux 


-^thresh 






[hr] 


[^] 


[xlO-^ m -V^] 


[TeV] 


1995/6 


Catanese et al. 1998 


20.5 


5.8 


1.7 ±0.5 


0.35 


20 Dec. 1995 


Catanese et al. 1998 


1.85 


5.3'' 


6.6 ± 1.9 


0.35 


1996/7 


Catanese et al. 1998 


24.9 


0.4 


< 0.82<= 


0.35 


Dec. 1997 


Aharonian ct al. 2000 


15.8 


NA 


< 0.29^^ 


1.0 


1997-2002 


Tluczykont ct al. 2003 


72.5 


4.4 


0.08 ±0.03 


0.8 


1998 


Konopclko ct al. 1999 


23.8 


3.3^ 


< 0.09'^ 


1.0 


2000 


Badran 2001 


3.1 


2.4 


1.1±0.P 


ftiO.4 



"^Statistical excess. 

''Part of the data listed in the above entry. 
^=99.9% C.L. upper limit 
•^99% C.L. upper limit 
'^Statistical error only. 
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Table 2: Details of the observations taken on 20 Dec. 1995. 



UTC^ 


El. 


Throughput^ 


7-rate'^ 


2:32 


64° 


0.77/0.78 ±0.08 


0.70 ±0.28 


3:34 


55° 


0.68/0.71 ±0.08 


1.04 ±0.37 


4:35 


47° 


0.52/0.57 ±0.05 


0.91 ±0.42 


5:36 


37°/36° 


0.45/0.41 ±0.04 


1.54 ±0.47 


Average 


1.14 ±0.20 



"Start time of the ON observation. The first three ON observations lasted for 28 minutes followed by an 

OFF observation. The last observation had a length of only 10 min. A complementary OFF observation, 
necessary for the spectral analysis, was chosen based upon similar observation conditions. 

''Relative cosmic-ray rate for ON/OFF observation, see text. 

^Gamma-ray rate per minute after Supercutsl995 (Catanese et al. 1998). 
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Table 3: Event statistics and flux in each energy bin for data set A. Upper limits are given 
at the 98% confldence level. 



Energy 


ON 


OFF 


ON-OFF 


S 


Flux 


[TeV] 


[events] 


[events] 


[events] 


W] 


[TeV-^m-^s-i] 


0.56 


63 


38 


25±10 


2.5 


(1.29±0.51)xl0-6 


1.12 


83 


63 


20±12 


1.7 


(1.27±0.77)xl0-^ 


2.24 


39 


42 


-3±9 


-0.3 


<3.91xl0-^ 


4.47 


22 


19 


3±6 


0.5 


<1.46xl0-^ 


8.91 


8 


7 


1±4 


0.3 


<3.62xl0-9 


Total 


220 


174 


44±19.8 


2.2 
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Table 4: Event statistics and flux in each energy bin for data set B. Upper limits are given 
at the 98% confidence level. 



Energy 


ON 


OFF 


ON-OFF 


S 


Flux 


[TeV] 


[events] 


[events] 


[events] 


W] 


[TeV-^m-^s-i] 


1.12 


55 


30 


25±9 


2.7 


(3.64±1.34)xl0-^ 


2.24 


86 


51 


35±12 


3.0 


(7.24±2.42)xl0-^ 


4.47 


35 


20 


15±7 


2.0 


(9.82±4.85)xl0-'^ 


8.91 


14 


6 


8±4 


1.8 


(2.11±1.19)xl0-9 


17.78 


7 


4 


3±3 


0.9 


<1. 19x10-9 


35.48 





2 


-2±1 


-1.4 


<1.43xl0-^° 


Total 


197 


113 


84±17.6 


4.8 





